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DIFFERENTIAL EQUATIONS, AND OF THEIR FINITE 
DIFFERENCE APPROXIMATIONS, WITH RESPECT TO THE 
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Abstract. In this article we study a controllability problem for parabolic and hyperbolic partial 
differential equations in which the control is the shape of the domain where the equation holds. 
Given a right-hand source term, the quantity to control is the trace of the solution into an open 
subdomain and at a given time. The mapping that associates this trace to the shape of the domain 
is non-linear. In this paper we first consider the continuous control problem and show an approximate 
controllability property for the linearized parabolic problem and an exact local controlability for the 
hyperbolic problem. Next we address the same questions in the context of a finite difference spatial 
semi-discretization. We prove a local controllability result for the parabolic problem, and an exact 
controllability for the hyperbolic one, applying the local surjectivity theorem together with a unique 
continuation property of the underlying adjoint discrete system. 
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Introduction. In this paper we consider parabolic and hyperbolic partial differ- 
ential equation on an open set C M", given a source term. This domain 51 fully 
contains a simply connected open set oj. We assume that the solution of the partial 
differential equation restricted to w at time T is accessible by means of measurement. 
The problem we address is a control and accessibility problem: Which solution can 
we reach at time T when the open set 51 varies? 

This paper follows the work of Chenais and Zuazua [7] about the controllability of 
an elliptic problem with respect to the shape of the domain. In our paper, the domain 
can vary in time and the equations are time-dependent, of parabolic and hyperbolic 
type. 

For the sake of clarity, we only study here the heat and wave equation with 
Dirichlet boundary conditions. Nevertheless, the method we use in this paper is 
valid for different parabolic or hyperbolic systems and directly applies for an elliptic 
equations. The same framework can also be used to study other boundary conditions, 
but the complexity of the problem increases. For instance for the same problem with 
Neumann boundary conditions require in the semi-discrete case a carefull study which 
is out of the scope of the present paper. 

The controllability problem is solved in the first section, and we show an approx- 
imate controllability property for the linearized parabolic problem and an exact local 
controlability for the hyperbolic problem. 

In the second section we address the same questions for the finite-difference space 
semi-discretized problem. Performing a finite-difference spatial discretization on the 
heat and wave equations, we address the same question as in the continuous-space 
case: is the solution of the semi-discrete equation controllable with respect to the shape 
of the (time-dependent) discrete domain?. The answear is yes, since we prove that 
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both the parabolic and the hyperbolic semi-discrete problems are exactly controllable. 

1. Control of the continuous wave and heat equations. In this section, we 
solve the continuous-time problem, and show that the parabolic linearized problem is 
approximately controllable, while the hyperbolic problem is exactly controllable. 

1.1. Position of the problem. Let £ be a parabolic or hyperbolic differential 
operator, / G i^(]R"), T > the time when we want to control the equation, and flo 
an open regular domain of R" . 

Let yo be the solution (when it is uniquely defined) of: 

Cyoix, t) = fix) V(x-, t)enoX R+ 

?/o(a;,t = 0) = Wx£W (1.1) 
yo{x,t)=0 V(a;,i) in aOo X R+. 

Let now V{jjq{-,T)) be a neighborhood of the function ?/oU 

(•,r) in H^iio). In 

this paper we are interested in finding small dynamical perturbations of the open set 
fio so that we can reach at time T any function yd G Hq{lu) n V{ya{T)). 

This problem is a control and accessibility problem, where the control is the shape 
of the domain. Many problems of this type are solved by optimization methods (see 
for instance [H El [21 [TU [161 [iTl [H]). In this case the problem is written in the form: 

i%\\yn-y4 

where U is the set of domains we take into account. Under suitable conditions on 
the set of admissible domains U, the existence of the minimizer can be guaranteed. 
Note however that this existence result in itself does not provide any information on 
whether j/^U = Vd does actually hold or not for a suitable choice of the domain fi, 
or even about the minimal distance between ynU ^-nd Vd- Therefore, optimization 
techniques will not solve the controllability problem under consideration. 

There are several versions of the controllability property that make sense in the 
situation under consideration and that we shall address. Indeed, the problem may be 
formulated at least in two different ways, as an exact or approximate controllability 
problem: 

(i) . The exact controllability: it is the case when one can find an open set n*{t) £ 

U such that 2/n-(r)|^ = ya- 

(ii) . The approximate controllability: for any fixed £ > 0, is it possible to find 

ri* (t) £ U such that for some norm on a certain set of functions we have 

II ynt{T)l~yd\\<el 

There is an extensive literature on exact and approximate controllability problems 
for partial differential equations (PDE) but very little has been done in the context 
in which the control is the shape of the domain. The interested reader may find an 
interesting introduction to these topics in J. L. Lions [Tl] and an updated survey 
article in E. Zuazua [201I2T] . 

1.2. Controllability of the heat equation. In this section, we consider the 
classical heat equation with a source term / e i^(K+,]R") on the open set whose 
frontier denoted by Fq is W^'°^: 
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dtya{t,x) — Ayo{t,x) = f{t,x) Vx e fto and t > 
yo{t = 0,x) =0 yxeilo (1.2) 

2/olanoW =0 yt>0 

This equation defines a unique solution ya , and we will perturb the domain flo to 
reach solutions in the neighborhood of yo{T)\^ in H^{lli) 

To this purpose, let us define the set U of admissible domains we take into account 
in the control problem. 

Definition 1.1. We denote by M^''-°°(R", M™) the set of functions k times dif- 
ferentiable o/ M" (in the sense of distributions), taking its values in R™, all the dif- 
ferentials being in L°°(]R", R™). 

We work in the standart setting for differentiation with respect to the domain (see 
for instance [l4 t fT6 t fT8 l fT7 t ll9j). The "admissible" perturbed open sets are the "small 
dynamical perturbations" of Oq. Mathematically speaking, let V a neighborhood of 
in the space of continuous and bounded functions L°°(R+, W^'°^{W\W)) f] C(]R+; , W^'°°{R"-,R")) 
such that for all (/3 £ V and for all t e [0, T], the maps {id + ^{t)) and {id + <f{t))~-^ 
are homemorphisms of H^^'°°(]R", M"). This set of (p functions is now denoted by W. 

For (f eV, let 

n^{t) = {id + ip{t)) {no) ^ {x + ip{t){x); X e no}. 

The perturbed equation is: 



dtyv{t,x) - Ayip{t,x) = f{t,x) V.T e n^p{t) and t > 
< y^{t = 0,x) =0 \/x eno (1.3) 

.y<p\T,it){t) =0 Vi>o 

We finally denote by A the non-linear map defined by: 



A 



W 



y'p{T) 



(1.4) 



and 7^(A) = {y^{T) \ip G W} the accessible states at time T of the perturbed equation. 

1.2.1. Existence and uniqueness of solutions for the perturbed prob- 
lem. In this section we use the variational formulation and show weak existence and 
uniqueness of the solution of the perturbed system. 

Theorem 1.2. The perturbed problem (|1.3p is equivalent to the following partial 
differential equation: 



\det{id + V^{t))\ 

y^{t = 0) = 

iV^ka = 



div{B{ip)Vy^) 



f o {id + ip) on no X R+ 

onno ^^-^^ 
Vt > 



where 



B{ip) ^ det{id + Vtp{t)) {V {id + ip{t))) 



(V(^d + ^(^)))* 



(1.6) 
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the star denotes the adjoint and 



{t,x) = y^(t,x + (p{t){x) 



Proof. We use the variational formulation of the equations l|1.3p and show that it is 

equivalent to the variational formulation of the problem (|1.5p . Let (j) £ ^R+, iJg (il^(t)) 

Let 4" the function defined on M+ x Qq by "^(t, x) ~ (j){t, x + ip(t){x)). 
We have: 



[dty^cf) ~ ^Vip'P) dtdx = / fcj) dtdx 



idty,p4> + yy^\74>) dtdx = / fcf) dtdx 

Using the change of variable x = r + if{r) with r G ^q. The determinant of the 
Jacobian matrix of the change of variable is det{id + Vip), so the change of variables 
theorem gives us the following equivalent expression: 

/ - (Vy^) o {id + (p) ■ (V(^) o {id + (p)} |det(id + V(p(t))| dtdx 

JS.+ X 0,0 

= f o {id+ ip)^\det{id + Vtp{t))\dtdx (1.7) 

JR+ xOo 

Let us now express (VG) o [id + (p) in function of v(Go {id + ip)y Let H = 
G o (id + Lp) and T = {id + ip). We have: 

dH _ d((GoT){y) 
dri dri 

^ do dxj 

■Hf dxj dvi 

= (VT)*VG 

where the star denotes the adjoint operator. Indeed, if x = Tr, 



/ dxi dxi 

dri dr2 

dx2 dx2 

\ dri dr2 



and hence 

dG 

VG=| ^ 

dx2 



Putting all these results together we get: 
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/ 

Jr 



R+xn, 



/]R+ xOo 

Finally, we get the following variational problem: 



Vy^ (V(id+V3(i)))* \/'i\Aei{id + Wip{t))\dtdx 
f o{id+ (^)^' |det(id + Vip{t))\ dtdx (1.8) 



/ {dtW^ \detiid+V^{t))\ - div(B(v?)) Vj?^} '^dtdx 

f o (id + ip)^ \det{id + Vipit))\dtdx (1.9) 



with 



B{ip)= det(id + V(^(i)) ( {V{id + ip{t))) 



{S/{^d + ^{t))y 



So the problem is equivalent to the partial differential equation on : 



^'^^ \dei{id + Vip{t))\ 

y^{t = 0) = 







Aw{B{Lp)Vy;^) ^ J o{id + ip) on VL^ x 

on f2o 
Vi > 



□ 

Theorem 1.3. The problem (|1.3p /las a unique solution in L'^{Q,T; Hq{^o)) D 
C([0,T],L2(17o)) 

The existence and uniqueness of solutions is an application of the Lions' theorem 
[21 theorem X.9] which we recall in theorem 1 1.41 

Theorem 1.4 (Lions). Let H := L^{no), V = H^{no) and V = H-\no) the 
dual space of V. Let f G L^((0, T), V') and uo G H . Consider the variational problem: 



{dtu, v) + a{t; u{t), v) 
u(0) 



{f{t),v) for almost all t e [0,T],Vu G H^i^o) 
uq on flo 

(1.10) 

Assume that a(t; u, v) satisfies the conditions: 

(i). the mapping t ^ a{t;u,v) is measurable 
(a). {u,v) I— > a{t;u,v) is bilinear for all t e [0,T] 
(Hi). \a{t;u,v)\ < M|lw||||w|| 
(iv). a{t;u,v) > a\\v\\l, - C\v\jj 

Then there exists a unique u such that: 

(I), u e L\Q,T-Hl{n^))nC{[Q,T],L\no)) 
(II). dtu&L^iO,T;V') 

(III), u is solution of the variational problem l|1.10p . 
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Proof, [of theorem II. 3| We consider the energy scalar product in defined by 



(/,5) = / f{x)g{x)\det{id + V^)\dx 



In the case of equation (|1.3p , we have 

a{t;u,v) = / B{ip)\7vW 



dx. 



/o„ v|det(^c^ + V¥>)L 

which is clearly measurable in t so the condition (i) of theorem 11.41 is satisfied. 

Furthermore we have B{ip) tends to identity when (/s — > in W so for small 
enough, | ||-B(</')|| - 1 | < 1/2. We also have |det(id+ V<^)| 1 when ^ in W so 
if we take ip small enough, ||det(i(i + V(/3)| — 1| < 1/2; so |det(id + V(p)| > 1/2 and 
||S(</')1I < 3/2 simultaneously. 

Eventually we get ,detH+v^)| ^o^"") 

We have: 



|a(t;u,w)| < 



V 



|det(id + Vv3)| 
|det(i(i+ Vv?)| 



dx 



<3/2||Vt;|| 
< 3||«|Ih.! IK^IIh.! 

Let now prove the last inequality on a of the assumptions of Lions' theorem: 



Jito V|det(id + Vv3)| 



1 

> - 



> ^ V 



|det(i(i+ V^)| 



So eventually theorem 11.41 applies and we have the existence and uniqueness of 
j/^, so existence and uniqueness of y^p in L^(0, T; ffg (Slo)) H C([0, T], i^(f2o)) (resp. 
L2(0, T; iioHf^^)) n C([0, T],L^{n^)). □ 

Remark Since the solutions of these equations exist, are unique continuous and 
defined at i = T, it makes sense to get interested in y^{T) and y^{T). 

1.2.2. Linearized problem :. In this section we define and study the linearized 
problem. The linearization is performed with respect to the parameter (p, and we will 
only consider the Lagrangian differential, which we denote y'^. 

Let us first recall the definition of the Lagrangian differential. 

Definition 1.5. The Lagrangian shape differential y'^ is the unique function 
(when it exists) satisfying 



VO (£ f7 y^lj-^ = yolsi + y'^\n + o{(p) 
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Remark We write A <£ B if and only if A CB. 

Proposition 1.6. The Lagrangian shape differential y'^ satisfies the following 
equation: 



dty'^it, x) - Ay'^it, x) = V(t, x) G R+ x fio 
y'^{t ^0,x)=0 Vxeno 

-ip.n ^ V(i, G R+ X dno 



(1.11) 



Proof. Let fl ^ il. We have in this open set the evolution equation 

dty^ -Ay^^f V(t, 6 R+ x f2 
Since yo satisfies the equation: 

dtyo - Ayo = / on IR+ x fj, 
we necessarily have the following evolution equation for y'^ on 11+ x il: 

dty'^ - Ay; = on R+ X 1^ 



(1.12) 



Now that we have the evolution equation of the Lagrangian shape differential, let 
us derive the boundary conditions. First of all, it is clear that 



y'{t = Q) = Q. 



Let us now compute y^p\^na■ Formally, we consider: 

yn{t)^pds 



for any ip G C°°(R"). We differentiate formally with respect to il, and we get the 
following formula (this calculation is rather classical, see for instance [U 6.28]): 



yn{t)4>ds 



n=no J do. 
y'^{t)ipds 

y'^{t)ipds 

y' {t)ipds 



ipUj.n — as 



(p{t).n{—i}j + yo — )ds 
on on 

ip(t).n——ds 
On 



because j/olasio = 0. Hence we have: 



(y'^it) + ^(t).n^)^ds = VV' G C°°(R2) 
So the following non homogeneous Dirichlet boundary condition follows: 
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and thus we get l|l.lip . □ 

1.2.3. Approximate controllability of the linearized problem. In this sec- 
tion we prove that the hnearized problem is approximately controllable using Holm- 
gren's theorem. It is clear that if c?A(0) was an isomorphism of iy^'°°(]R", R") on 
-ff^(R"), then we would have using the local inversion theorem a property of exact 
local controllability. The same property would be valid if dA(0) was surjective. But 
the regularity introduced by the parabolic form of the equation impHes that we will 
never have this controllability property. Indeed, the regularity of the solution is linked 
with the regularity of the source function /. The most interesting property we can 
expect is the density of A/'n W in N, or for the linearized operator c?A(0). 
As a consequence, in the sequel, we try to prove a property of approximated control- 
lability. We prove that the function dK(Q) actually has a dense image, which is an 
interesting property by itself. But nevertheless, even with this property we are not 
able to prove that A has (even locally) a dense image. 

So this method fails to prove the approximated controllability. This is the reason why 
we address in section [2?T] the same questions for the space semi-discretization. 

Holmgren's theorem [8] is based on the propagation of zeros and uses the charac- 
teristic surfaces of the underlying equation. 

Lemma 1.7. The characteristic surfaces of the heat equation are the hyperplanes 
t = constant. 

Proof. The heat differential operator we consider is 



So the associated characteristic polynomial reads: 



k=l 



P{T,X) ^T-\X 



|2 



and its principal part is P2{T,X) — \X\'^ — 0. The solution of this later equation is 
\/i = 1 . . . Xi = 0, and its direction (1, 0, . . . , 0), so eventually the characteristic 
surfaces are the hyperplanes T = constant . □ 

Proposition 1.8. Let flo an open subset o/R", whose boundary Tq is regular 
(say , or difjerentiable with bounded differential). Let 7 a non-empty open subset 
of To. 

Then any solution of the equations: 

'-dfu{t,x) - Au{t,x) = yt>o,xeno 

u(i = 0,a;) = Vx e flo , , 

(1 13) 

u\raix,t) =0 Vi > 0,x G To ^ ' ^ 

i^l,=o Vt>0,xe7 

is null on flo 

Proof We use Holmgren's theorem to prove this proposition. Let P{D) the heat 
differential operator : P{D) = dt — A. Lemma fTTZl ensures us that the characteristic 
surfaces of P{D) are the hyperplanes t = constant. 
Furthermore, we assumed that : 
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{P{D)u^O V(t,a;) G R+ X r^o, 

■"lrxiR+ =0 on To X K,+ , 
at; 
dn 



I 1^ = V(<,x)eM+X7 



First of all, let us write the variational formulation satisfied by the function u € 
C{0,T;H^{no)) solution of ifTT^ll : 

yveH\no),f {dtuv + VuVv)dx~ [ ^vds^O. (1.14) 
Jno Jto on 

Let now u be a continuation of u on an open subset close to flo defined by 

= f7o u 

where ^^(7) is a small perturbation corresponding to a smooth modification of the 
boundary of on 7. Let u G i/Q(17o) be the function defined by : 



u = u on ilo 
-u = on 0(7). 



(1.15) 



We have the following relations: 



9tw|n(7) = 
VM|a(^) = 
,§| =0onaf^(7) 



and hence we have 



r f du 

yveH\n{-/)), dtuv + VuVv- —vds^O. (1.16) 

Jn{-f} Jdni-(} 

Now, using l|1.14p and (|1.16p . we obtain that u satisfies the following variational 
problem: 

yveH^no), [dtuv + VuVv- [ _^vds = 0. (1.17) 

At this point, we are in the field of application of Holmgren's theorem, except 
that the subsets and ^^(7) are not convex. Nevertheless, we can find a partition of 
flo in open discs. Those discs are convex, and the intersection between two discs is 
also convex so we can apply Holmgren's theorem on those discs and then the theorem 
is proved. 

More precisely, let B be a disc from^is partition. Assume that B = BiU B2 and 
u\bi =0 (this is always possible when fio is connected). 

• Ai = M.*^ X Bi and fl2 = IR+ x B2 are two open subsets of R"+^, convexes, 
and such that £ ^^i 

• P{D) = 9t — A is a constant coefficients differential operator. 

• each characteristic plane P{D) intersecting U2 is an hyperplane t — constant 
which intersects always also wi. 
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• u solution of P{D)u — and satisfies also = 
So Holmgren's theorem applies and we deduce that u = on f^o, in particular on fio 
on which we have u = u, so we have indeed u = on fio □ 

Thanks to these results, we can prove the density result previously announced for 
the continuous problem: 

Theorem 1.9. Assume that there exist a non-empty subset o/Tq on which 

VtG[0,T], ^^0 

Then TZ = {y'^{T);ipe W} is dense in L'^iflo). 

Proof. We prove that TZ-^ ~ 0. To this purpose let us consider g G i^(rio) such 

that 



guj ~ 0, \/uj ^ TZ 



On 



and let us prove that necessarily g ~ 0. By definition, TZ = {y'^{T)\i~p G W} so the 
later condition is equivalent to : 



"^^eW gy'^{T) = (1.18) 
Let us define the adjoint state (t> G H^{Vl) associated to g: 



-9(0 — A0 ~ g ® 5t=T 

(j){t = 0) =0 

, iro =0 
with by definition 

yv eCiO,T;H\no)){g®St=T,v) = [ gv[T)dx. 

Thus we have to show that = i.e, g — 0. We have: 

dni) =^ G W{-M - A0, y;) = / gy'^iT) = 

J do 

We now integrate by parts and get 

+ I V0V,; - / 



f.n— — —atas 
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= 5 (St=T 


<P{t - 0) 


= 




= 




= 



We deduce that for all t in [0, T], we have = on Fq. 

We assumed that on the non-negligible subset 7 £ Fg we had 7^ V< G [0, T]. 
Hence the adjoint state 4> satisfies the following equation: 



(1.19) 



From Holmgren's theorem [8l theorem 5.3.3], the zeros propagate instantaneously 
on all the zone in contact with the characteristic surface crossing 7, so the solution (j) 
is zero for all t G R"*". Indeed, 

• for all t < T, the equation satisfied by <j> is the homogeneous heat equation 

9t0 - = 

so it is from the propertv ll.81 

• for all t > T, the equation satisfied by is also the homogeneous heat equation 

so it is null by the same argument. 
Lastly, the condition 

—dt4> ~ Acf) = g ® 5t=T 
is equivalent to say that the jump of 4> at T, 

[|0|] = 0(T+)-0(r-) 

is equal to g. Indeed, assume that g is piecewise continuous with a discontinuity at 
t — T. We write the variational problem associated to (j). Let v S C(5"(]R"'", i?o (fJo)) 

{dt(j)v + \7(j)Vv}dxdt = / gv{T)dx 

R+xfio Jn.0 

I [ (pdtv + [ [ (f>dtv + [ V<j)Vvdxdt = [ gv[T)dx 
Jo Jua -It Jn„ Jn+xna Jna 

{4>{T+) ^ <ly{T-))v{T)dx 



{<i)dtV + \/(l)\'v}dxdt^ / [g - {(l){T+) - (j){T-))]vlT)dx 
R+ X no Jn.0 

So g is equal at the jump of 4>. Here, 4> is constant, equal at so necessarily <? = 0. 
(QED) □ 

1.2.4. Conclusion. In this section we have proved that the linearized problem 
is approximately controllable, i.e. that it has a dense range on i7g(rio). This result 
does not allow us to go any further. 

Consequently, even if the density of the range of the linearized operator might 
be of independent interest, it does not have any relevance, a priori, for solving the 
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original problem in which the nonlinear mapping is involved. It is precisely for this 
reason that we will turn our attention to the discretized problem in section 12.11 

Nevertheless it was clear since the beginning that the exact controllability would 
not hold, because of the regularity of the parabolic operator. This is not the case for 
the hyperbolic operator, and we will show that an exact controllability property will 
hold. 

1.3. Controllability of the wave equation. In this section we study the con- 
trollability of the wave equation with respect to the shape of the domain. We prove 
here that the hyperbolic problem is exactly controllable. 

First of all let us define the set of admissible domains U for this problem. We are 
quite free on the space to chose for these perturbations, the only requirement is to 
have sufficient diversity so that the problem will be controllable. We do not look for 
minimal conditions, but just for sufficient conditions on the perturbations we admit. 

First of all, we require that the admissible open sets should satisfy the same 
regularity properties as the heat equation: 

W = L°°(Il+,V)nC(]R+,V). 

Moreover, since the "information" has a finite speed in the wave equation, the 
perturbations we define should be slow enough to ensure that the problem is well 
posed (the perturbation should stay in the characteristic "dependence cone" of slope 
1). More precisely, the time differential should be in L°° with norm less than 1. The 
perturbation (p should be in W, be time-differentiable and in the sense of distribution 
and its differential should be in L°° and ||9t(/3||L^ < 1. 

On the other hand, at a technical level, it is easier for us to assume that the 
perturbation is C^, so that the existence and uniqueness for the perturbed problem 
is ensured. Finally, to be able to define the space where we will invert the functions, 
we assume that the transformations are uniformly bounded by a given a. 
Eventually, we define the following set of admissible transformations: 

X = {(^ e Ci(]R+,V)nL°°(R+,V); ll^lloo < a, dt^ ^ {M + ,V) and ||9t(^||oo < 1 }• 
For e A^, we denote: 

n^{t) = {id + ip{t)) (r^o) ^{x + </?(t)(x); X e r^o}. 

Let T = 2(diam(f^o) + a) and B{y, R) ^ {x e R"; \x - y\ < R} la ball of center 
y and of radius R. Let also 

= {dna + B{Q,a)). 
= r^o + B(0,a). 

This is the "reachable neighborhood" of Tq at time T in R". When we will take a 
perturbation in the space we consider, all the perturbed open subsets will be inside 
the open set . Let us finally denote H^^^Qo) = {i?o(T);5T C X^} and ^'^(r^o) = 
{L%T);dTcX^} 
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Let / e L^{K+,H^\W')), y" e H^i^o) and e ^^(f^o) we define as the 
solution of the equations: 

'd^y^it, x) ~ l^y^{t, X) = /(t, x) > 0, a; G ^^{t) 

y^{t = 0, x) = y^{x) Va; G fio 

9t2/^(< = 0,x) = y^(x) Vx e fio 
,y<^(t,a;) = > 0, x e V^[t) 

The question we address is to characterize the set of traces at i = T of the 
solutions of this problem when is an admissible transformation: 

n{M) = {{y^{T),dty^{T))-^>^ M\. 

Let also j/o be the solution associated to the trivial perturbation </? = 0. 

Let us finally denote: 

■ \^ ^ {y^{T\dty^{J)) 
and 7^(A) = {{y^{T),dty^{T)) G X}. 

1.3.1. Existence and uniqueness of solutions. The same method as the one 
we used for the heat equation yields to the following equation: 



/ 



{5f2y^|det(id+ V(^(t))| - div(B(</j))Vy^}M(i2; 



+ xOo 



! o{id^ if)^\A(i\.{id + \I^}(t))\dtdx (1.22) 



where 



= det(id + V^(t)) (V(i(i + (^(t) )) 



V(id + </?(t)) 



We can see that -8(1^9) is symmetrical. It is positive because it is a perturbation 
of identity (when ^) is small enough). We now work on with the dot product : 



((a,6), (c,(i))Hi(Oo)xL2(Oo) = / B(<<5)Va.Vc+ / hd\As±[id + SI ip{t))\dx . 

Here again we use Lions' theorem 11.41 to prove the existence and uniqueness of 
solution for the perturbed system. 

Indeed, with Lions' notations we have: 

a{t;u,v) = (u(t),u)^i(no). 

Thus a is bihnear, symmetrical, and continuous. We only have to prove that it is 
. The time dependence is directly linked with the function Lp. Taking t ip(t) in 
C^, we are sure that t ^ a{t) is . So under this condition on we have the first 
assumption of Lions' theorem. The other hypothesis is the same as the one we used 
for the heat equation and the same proof appHes. 

Hence Lions' theorem [l .41 ensures us that there exists a unique y^, so a unique y^p. 
The solutions exist until t ~T and are continuous so it is possible to be interested in 
y^{T) and y^{T). 
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1.3.2. Linearized problem. In this section we compute the hnearized problem 
using the Lagrangian shape differential we defined in the heat equation section. 

Proposition 1.10. The Lagrangian shape differential satisfies the following 
equation: 




(1.23) 



1.3.3. Exact controllability of the wave equation. In this section we prove 
that (iA(0) is surjective from M to i/o^(ilo) x L'^(rio). The local surjectivity theorem 
will give us as an application the exact controllability of the problem 

Theorem 1.11. The linearized function dA{0) is surjective from C'^{[0,T],Wh) 
onto i?o(^o) X L^irio) 

Proof. The proof of this theorem is based on classical results on the zero con- 
trollability of the wave equation, proved by Lions in \W\ [T2] . In this paper, Lions 
considers an open bounded subset fl of R" with a smooth boundary F and a con- 
trol time T > 2diam(ri), 7 C F satisfying some geometrical control conditions. The 
problem was to find a function v such that the solution of the equation 



ifu - Au 





on X 


(o,r) 


uiO) 


= Mo 


on Q. 




dtu{0) 


= «1 


on Q. 




u{x, t) 


= v{x,t) 


on 7 X 


(0,T) 



was such that yv{x, T) = dtyv{x, T) = 0. 

Let us return to the proof of our theorem, 
defined by: 



We recall that the function dA(0) is 



dA(0) : 



^ ^{y'^(T),d,y'^{T)) 



where the function y'^ is solution of : 

dy^~Ay'^ = 

y;(i = o) = 

dty'^{t = Q) = 

Let us assume that ^ > e > on a subset 7 C Fq satisfying the geometric 
control condition of Lions [lO]. In this case, —ip.n^^ can take any value in ^^(7) and 
we can exactly control this equation for all time T > 2( diam {Qq) + 2a). 

More precisely, there exists a certain ^ € M vanishing outside 7 and such that for 
all pair {y^ ^y^) of initial condition, we have yxp{T) = y'^ and dty4,{T) = y^. Indeed, 
there exists ip vanishing on the boundary of flo outside 7, such that the problem 
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satisfies y^(T) = and dty^(T) = 0. But since the wave equation is reversible, 
changing t into t = T — t, we obtain that y^ir) := Y{t) satisfies : 





AY 







T) 


- y° 


Yit^ 


0) 


= 


dtY{t 


= T) 


= 


dtY{t 


= 0) 


= 






= -v^ 



So Y is the solution we are searching for. 
□ 



1.3.4. Control of the wave equation. In this section we prove a theorem of 
exact controllability of the wave equation with respect to the shape of the domain. 
Theorem 1.12. Letya the solution of the unperturbed problem: 

dfyo - Ayo = / on flo 

yo{t = 0) = on VLq 

dtyoit^O) = y^ onVto 

. yolrowW = Vi>0. 

There exists a neighborhood of {yo{T), dtyo{T)) in Hq{Hq) x i^(r2Q ), denoted Af such 
that for all A G Af, there exists (p <E AA such that A = A((/3). 

Proof. We proved that A was differentiable and that its differential at is surjec- 
tive. So the local surjectivity theorem (see e.g. Luenberger [l^ )proves theorem [1.121 
□ 

2. Control of the semi-discrete wave and heat equations. In the section 
[H we addressed the problem of the controllability of the heat and wave equations 
with respect to the shape of the domain and proved section that the wave equation is 
locally exactly controllable with respect to the shape of the domain, that the paraboHc 
problem was not exactly controllable and that its linearization was approximately 
controllable. 

For this reason that we shall turn our attention to the discretized problem. There 
are two good reasons for doing this. First, it is relevant from a computational point of 
view and, from the mathematical point of view, since the range is finite dimensional, 
the density of the range of the linearized operator will imply its surjectivity and 
therefore allows to apply the inverse function theorem. 

In this way, in the context of the discretized finite-dimensional equation, we shall 
prove a local exact controllability result. 
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2.1. Semi-discrete heat equation in a square. In this section we consider 
the semi-discrete heat equation on a square [0, a] x [0, b] G M^, with a/b e Q (so that 
we can use the mesh size in the two directions). We discretize this set with a step h. 
The domain we consider is now 

a, = {m={ih,jh);{i,j) e {0,...,A/} X {0,...,7V}}. 

Definition 2.1. Let m = {ih,jh) e (Z/i)^ and the discrete neighborhood of m 
defined: 

B{m) = {{kh,lh)- {kj) = (z, j), (t - 1, j), {i + l,j), J - 1), + 1)} 

The set of strict neighbors of m is B{m) ~ B{m) \ {m}. 

Definition 2.2. The discrete interior of flh is defined by as: 

the discrete boundary of flh by 

F/i = nii\ Qh, 

and the exterior of flh as: 

Remark Those three sets form a partition of {"LhY . 

We assume that the only free part of the boundary is = 0}. This means 

that the only moving part of this set is the only part of the boundary that can move 
to control the solution. Furthermore, the free points of the boundary will move only 
along the normal to this boundary. 

Definition 2.3 (Function spaces). We denote T{X) the set of maps X ^ 'B. 
and J-i){X) the set of maps in J-iX) vanishing on the boundary of X. 

We consider in this paper time dependent maps, taking values in ^{flh). In 
particular we will use C([0, T], jr(r2/j)), the set of continuous functions [0, T] ^ J-'{flh) 
and the spaces U'{Q,T;T{Vl.h))- 
Remark The set will be identified to R*^^. 

2.1.1. Unperturbed state. DEFINITION 2.4. Let A be the finite difference 
operator with Dirichlet boundary conditions defined by: 

(j) — > A(j) 
where 
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The reference state u we consider is the solution of the equation 

_ , , , f dfU + Au = F „N 
u G >o(' 'ft ) such that < „s (2.2 

^ ' [ u[t = 0) = Mo 

with F G J-{nh) depending on the source term / of the continuous initial problem. 
If / is continuous, then F will be: 

Vm eflh , Frn = f{m). 

If / is not continuous, for instance is / € or H^, then we can change the value of 
/(to) by a mean value of / on a neighborhood of to. 

2.1.2. The perturbed problem. As in the continuous case, we are interested 
in small perturbations of the shape of the domain Qh- Changing the shape of flh 
consists in moving continuously the nodes of the mesh corresponding to a; = 0. On 
this new subset, the finite difference Laplace operator is modified as follows: 

Definition 2.5. We consider the set {Vj-J = l..iV-l} of vector fields ilh 

by: 



y? e {i..iv-i}, 



yj(TO)-(0,0) ifm^{0,jh) 
F,(to) = (1,0) ifm^{0,jh) 



Let Wh be the vector space spanned by the family (Vj)jg{i.. jv_i} . The perturbations 
we consider in this problem are in the set: 

N-l 

Wft = {JZ ^-^jW^^;* ^ ^^(0 e L°°{]R+;Wh) n CCR+iWh) 

and such that sup || Aj ||oo < 1/2} (2.3) 

j = l..N-l 



Remark Note that the perturbation has the same magnitude as h. This does not 
allow us even in the better cases to have the continuous case as a limit case since the 
perturbation tends to the trivial condition as the mesh becomes finer. 
Definition 2.6. Let us define the first layer of interior nodes: 

r^ = {(i,j) ;je{i,..iv-i}}. 

Let us now define the perturbed Dirichlet operator. 
Definition 2.7. Let ip{t) = J^^S^^ >^3{t)hVj G Wh- 

The operator A{ip) : To{n,f^) J-'i^.h) is defined as: 

r^[40(TO)- E cj^ip)] VTOGfi/Ar^ 

l7r?P(l + T+I-w)'?^(ij) - 2+^'?^(2j) - -'^(ij-i)] /or TO = (l,j) G ri 

(2.4) 

Proposition 2.8. The operator A{(p) is bounded for all (p G Wh- 
Proof. 
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For merih \r\ we have [A{ip)(j)]m = [A0],„ M'^(t>im)-J2p(zB{ra},p^m 
so we have: 



<^ll^ll 



■ E I 

B (m) .p^m 



For m = G r^, we have 

11 2 



/l2^ ^ ' 1 + Aj(t) 

and hence we have: 



Since we are in J-{flii) which is a finite dimension vector space, all the norms are 
equivalent so the operator A{ip) is bounded on V. □ 



The perturbed state 

Definition 2.9. With the definitions introduced in the latest section, we define 

the perturbed state u^{x,t) g ^{Qh) as the unique solution in C{[0,T], T{Qh)) of the 
semi-discrete problem: 

idt^ + A{(p)^ ^^^^ 
yu^{t = 0) 

This solution exists and is unique. 

Proof. Here again we write the problem as: 

idt^ = -A{lp)^ + F ^^^^ 

|u^(t = 0) =ut) 

So it is a Cauchy problem 9(1? = X(lE,t), with X^Xji) = —A{ip{t))'x + F{t), 

in R*^^"^) ~ The function XCliCjt) is measurable with respect to i; it is 

Lipschitz in ~x from proposition [231 So we have local existence and uniqueness of the 
perturbed state, which will be continuous in its definition domain. □ 
Proposition 2.10. The perturbed state is defined for all t > 
Proof Indeed, the perturbed state is solution of the equation 

dtu^ ^ F - A{ip)u^ 
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Let us use Gronwall's lemma to prove that the function cannot blow up in finite 
time. 

Wdt^W < \\A{^K\\ + \\F\\ 

< a\\u^\\ + /3 

If we assume for instance that F e L°^{[0,T],Wh) or F e C{[0,T],Wh), uniformly 
bounded by /3 on [0,T]. The a-constant used in the later equation is the bound of 
the operator A coming from the proposition 12.81 
Let p be the solution of the equation 

= ap + (3. 

= K(i-o)|| 




a) a 



The solution of this equation is 

p(t) = (p(0) + ^)e° 
So Gronwall's lemma ensures us that 

KWII<(lk^(0)ll + ^)e"*-^ Vt>0 

Eventually, the function u^it) has a finite upper bound for all t > so it does not 
blow up in finite time, so it is definite for all positive finite time. □ 

2.1.3. Controllability of the semi-discrete heat equation. As in the con- 
tinuous case, we consider the map 

Where u^{T) is solution of the equation l|2.5p . Let Zd = u{T), where u is the refer- 
ence state (|2.2p . The problem we address is to find a neighborhood V(0) S V of the 

reference domain and V{Zd) G T{^h) of the trace at t = T of the reference solution 
such that V{Zd) C Ah(V(0)). 



In the finite-dimension spaces where the problem is now set, we prove in this 
section that this is possible. To this purpose, we will use the local inversion theorem. 
First of all we will prove that Ah is differentiable in the neighborhood of the origin, 
and that dAh{0) is surjective. Then we will use the adjoint state technique to prove 
that the surjectivity of dAh{0) is equivalent to a pool of conditions the semi-discrete 
adjoint should state satisfies. Finally, the we will prove the controllability property 
proving those conditions on the adjoint state, which happen to be a property of 
discrete unique continuation. 

(i). Differentiability 



Let us denote tr the trace operator: 
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Then A/j is the composition of the trace application and the map U : (p — > u^. 

The trace function is hnear and continuous. So we only need to prove that U is 
Frechet-diflerentiable in 0. 

Proposition 2.11. V(p g Wh , ^ip e Wh, the map ip is differentiable in 

if in the direction of tjj, and the Gateaux differential {DGAfi{Lp),ip) , denoted v^{ip) is 
solution of the differential equation: 




Proof First let Wx 



Uip+X'4> 



X 



{2.. 



The function Da^h is the limit, when it 



exists, of Wx when A tends to 0. 

(i). Necessary condition: We assume that this limit exists, denote it v, and we 
compute the equation this limit satisfies. The function u^+xtp satisfies the 
equations: 

\dtU^+X4, + ^(<<5 + AV')w<^+A^ ^ F{t) 
\u^+X4,{t = 0) =u^ 

We differentiate this equation with respect to A at A = 0, and we obtain the 
equation v, 

jdtv+ {A^,'ilj)v^ + A{ip)~u =0 
[^(t==0) =0 

So we deduce that if the diflterential of Wx exists, then it is the solution v^p{'ili) 
of the ordinary differential equation (|2.8p : 

] dtv^{ilj) + A{(p)v^{il}) = -{A^,ijj)u^ 



= 



(ii). 



Sufficient condition: We show that the solution of the ODE l|2.8p is the limit of 
Wx- Indeed, the map {Xj)j=i,,,N-i i— > A{(j)) is C°°, so the differential {A^,ip) 
is defined, and furthermore {A^,ijj)u^ is L^(0,T). So Cauchy-Lipschitz gives 
us the existence and uniqueness of the solution, and we can prove that it is 
defined for all time using Gronwall's lemma. 

Let us now show that the function ~v'{t) defined is indeed the limit of Wx 
when A 0. Let Wx ~ v :~ £\. We have: 



ex = dtex + A{ip + \ij)Wx 
A{ip + Xip) - A{ip) _ 
A 



- A{ip)v 



Since ^ e in ip, the right hand of the equality is ©(A), as well as 
{A{(p + Xip) — A{lp))v which comes from the left hand of the equality. So 
eventually, ex satisfies the equation: 
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So we get that 



From Gronwall's lemma, we have ex 
Indeed, let p be the solution of 



dtSxW < a||eA|| +A/A 
when A - 



0. 



P 

Pit - 0) 



ap + \M 




so we have: 



dtipe "*) = (p - ap)e~ 



p{t)e- 



P(0) 
Pit) 



A 



le 

AM, 

a 
M 



-(e-"*-l) 



1) 



And eventually, using Gronwall's lemma, we have 

'M 



||£AWII<A(-(e"*-l)) 
^A(^(.--l,), 



estimation which is independent of t. So e\ converges uniformlj0 to when 
A ^ 0. So the limit of W\ exists and is indeed the function v solution of 
(ESI). 



□ 

Lemma 2.12. Let f : R" R" such that f is Gateaux- differentiable in all 
directions in the neighborhood of a point xq, with continuous differentials. 
Then f is Frechet- differentiable at xq and 



This lemma is quite famous in differential calculus, and can be generalized in 
infinite dimensional cases: 

Lemma 2.13 (Generalization). Let N : X ^ Y , and 

DgN : X CiX, Y) 

the Gateaux- differential of N . If DgN is definite and continuous in the neighborhood 
ofO, then N is Frechet- differentiable at and DgN = DN, i.e : 

\\N{x + y)-N{x)~DGN{x)y\\ 

hm — = U 

\\y\\^0;yex \\y\\ 

Proposition 2.14. The differential of with respect to 4>, 

V — * DGAhiip) 

^We could have used the Cauchy-Lipschitz theorem with parameters to prove this last point 
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is continuous. 

Proof. DQAfi{(p) : tp — > 'v where v is solution of: 

V(i = 0) = 
Hence "u* is solution of 

5f = F^^^it, ~v) 

with 

F^^^{t,lr) = -A^~x - {A^,'\l))y^. 
We note that the function 

is continuous, since Lp is continuous. Cauchy-Lipschitz' theorem with parame- 

ters gives us that i-p ^ y^'vi also continuous. 

Eventually, the differential {A^,%1}) is continuous in ip: indeed, A^ has a rational 
variation in (p, and Lp are such that those rational fractions have no singular point, so 
the dependence in Lp remains continuous. 

So Cauchy-Lipschitz' theorem with parameters gives us the continuity of in kp. 

□ 

Theorem 2.15. A/i is Frechet-differentiable atO. 
Proof. We already proved that: 

• the Gateaux differentials in all directions of V exist (PropETTT]) 

• those differential are continuous (Prop. I2.14p 

So lemma [2TT31 gives us the Frechet-differentiability of A/i at 0. □ 

(ii). Adjoint state technique 

We recall that ip € 1^(0, T; Wh). We have: 

A/i, defined in (|2.7p is the compound of the trace function 

I F F{T) 

and the map denoted U defined by: 

^ . \ C{[Q,T]-Wh) — C([0,r];.F(4)) 
I (p — > 



where u^p is solution of the equations: 



dtu^ + A{ip)u^ = F 
u^{t = 0) ==0 
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Furthermore, recall that 
and yip is solution of 

r dty^ + A{{))y^ = -{A{Q),ip)uo 
\ y^{t = 0) = 

Eventually, we denote Y : ip — > y^ (so we have dKh{Q) = tr\t=T ° Y)- 

The aim of this section is to show that the differential of at is surjective. 

Let c S T{Q.h). It is clear that the adjoint of the trace map, tr*c, is c6t=T- We 
now prove that the map dAft (O) is surjective. To this purpose, we use the adjoint 
state technique. 

dAh{0) is surjective ^ {c £ T{dh);y^ £ Wh{dAh{0)^,c) = 0} = {0} 
We also have: 

^ipeWh {dAh{0)ip,c) = 
^ yipeWh {tr\t^TiY{ip)),c) =0 



^ G Wh{Y{^), trll^ric)) = (2.9) 
Furthermore we have Y{ip) = y^ is solution of the differential equation: 

dty^ + A{Q)y^ = -A'^yo 

Definition 2.16. The adjoint state associated to and tr\^^ j,{c) is the unique 
solution X of the equations: 



Remark 



~dtX + AX =tr\l^j,{c) 
X{t^O) =0 



{dty^ + Ay^, T) = {dty^.T) + {Ay^, T) 

= -{y^,dtT) + {y^,A*T) 
= {y^,~dtT + AT) 



(2.10) 



(we used that A is self-adjoint {A* ^ A) ) 

Using this definition we replace in (|2.9p tr\^^rp{c) by its expression in function of 
the adjoint state X defined in l|2.10p . 



Hence we have: 



yip€Wh{Yiip),tr\t=A<^))=0 

^ e Wh{Y{ip), ^dtX + AX) = 

^ Vcp e Wh {dty^ + Ay^,X) = Q 

^y^eWh{A'^yQ,X) = Q 
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So the following theorem holds: 

Theorem 2.17. The differential dAfi{Q) of at (p ~ is surjective if and only 
if we have the following unicity property: If c £ J-{Vlh) is such that 



yo) o = 0, Vv3 e W,, (2.11) 



where X is solution of : 

~dtX + AX = tr\;^j.{c) 
X{t = 0) = 

Then necessarily c — 0. 

(iii) . Calculation of the differential of ^ at 

Proposition 2.18. Letj e {l...iV-l} andipe T{dh). For all n € C([0,r],R), 
we denote {A'q, ^{t)Vj) the differential of A at in the direction fi{t)Vj. We have: 

[{<,Kt)vMn. = l%,,,^ ^, . (2.12) 

{2^(2,1) - 2'/'(ij)) if m= 

Proof For all j G {l...iV — 1}, the point (0, j) of the boundary has a unique 

neighbor in flh, which is (1, j). From the definition 12.71 of A{ip), for all /i : t ^ /i(t) 
C([0,r];]R), we have : 

Vm en,, , m ^ (1, j) , [A^(t)v^.0]™ = 

and we have also 

1/1 2 \ 

[MKt)V])<P]{h3) = (^2(1 + - 2 + ^(t) '^(^'J) " "^(iJ+i) " '^(iJ-i) J ■ 

Hence we get 

[A{^,{t)vMii,) - [Mh,n = ^[2(1 + - 2)0(1.,.) - „ 

= 7^ [2(1 + 1 - Ai(t) + oi^i) - 2)0(i,j) - (1 - - 1 + o(/^))</>(2,,)] 

□ 

(iv) . The condition X\-pi =0 

Definition 2.19. The function Y <= C{[0,T],T{yLh)) satisfies the discrete non- 
degeneracy condition if and only if: 

v< > 0, V? e {1, ...,N - 1}, iy(2,,) - 2^(1,,) ^ 0. 

Remark This condition can be seen as an finite difference approximation of the 
condition 7^ 0. Indeed, let y G C^(V(x)) where V{x) is a neighborhood of x. 
Assume that y satisfies y{x) = 0. Then performing a Taylor expansion, we get: 

y'{x) = ^[2y{x + h) - ]^y{x + 2h)] + o{h). 
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So the condition 2^(2 j) ~ is the finite difference expression of the continuous 

condition 7^ we assumed in the continuous case, in theorem ll.91 
Note also that in the continuous case, we only need to assume that || does not 
vanish on an open set of the boundary, and not all along the boundary. Here we need 
to assume the non-degeneracy condition all along the boundary of flh to prove the 
discrete unique continuation. 

Proposition 2.20. Let us assume that the reference state satisfies the non- 
degeneracy condition all along the boundary offlh- Then we have: 

X\ri = 0. 



Proof. The relation (|2.1ip gives us : 
Proposition [2118] gives us therefore that Vj £ {1, .., iV — 1}, 



. ym erth , m ^ (1, j) , [A'„{^l{t)VJ)(P]„^ = o 

so we eventually have 



Thus 



1^ J2 KtMoiKt)Vj)yaU{t)X^{t)dt ^ A^(t)^[i2/o(2,j)-2yo(l,j)]X(i,,)(t)di = 

So we deduce that for all t > 0, ■^[yo{2,j) - 2yo(l, W = 

Since ^yo{2,j) — 2yo(l, j) never vanishes, we have: 

= Vt > and Vj e {1, TV - 1}. 

□ 

(v). Unique discrete continuation 

The aim of this section is to prove the uniqueness condition l|2.1ip . From proposi- 
tion [5201 and under the discrete non-degeneracy condition on the reference state yo, 
we have X\y-i = 0. 

We now prove this condition using the adjoint state defined by (|2.16p : 

~dtX + AX = tr\;^^{c) 
X{t = 0) = 

implies that X identically null, so trl^rp[c) = and c = 0. To this end, we study the 
propagation of the zeros of X on from its boundary, as we did in the continuous 
case using Holmgren's theorem. The main difference is that the propagation of zeros 
in the continuous case is a global property, whereas it is a local property in the discrete 
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case. This is why we need to assume the non-degeneracy condition all along the the 
boundary of flh- 

Theorem 2.21. The unique solution of the equations 

' -dtX + AX = tr\;^j,{c) 

X\ri = 

is X = 0, and so c = 0. 
Proof. 

(i). First we are interested in the equation ~dtX + AX = tr\^^rp(c). On the sets 
[0,T[ and C!o[, the equation simply reads —dtX + AX = so we have 
existence, uniqueness and continuity of the solution X in these domains. The 
right hand term can be interpreted as an imposed jump condition at time 
t = T. Indeed, let us write the variational formulation of the problem l| 2. 13p : 

Let V € C^(R"'", .F(r2ft)). The variational formulation reads: 

/•oo 

/ -dtXrnVrn + {AX)m..Vmdt = Vm{T)Cm 

i.e : 

/•oo 

/ y2 dtXmVm + {AX)„,.Vmdt = Vm{T)Cm 

The imposed jump at t = T reads: 

-dtX + AX = Q 
X{t = 0) = 

X|ro = (2.14) 

where we denoted [\X\]{T) = X{T+) - X(T-) the jump of X at i = T. The 
variational formulation of the problem reads: Vu S C^(R+, 



/•oo 

= / -dtXmVm + AXm,.Vmdt 

JO 

pT poo 
= / Y -dtX^V^+ AXrn-Vmdt+ / ^ "dtX^nVrn + AX^-V. 

j-T POO 

= / V] X^dtVm + AXjn-V,ndt + / y2 XmdtVm AXm-Vmdt 
Jo „ JT 



- J2 {Xrn{T+) - X^{T-))Um{T) + Xm{0)v^iO) 
pec 

/ V X„AVm + AXrn-Vmdt - ^ [l^m |] (T)«„, (T) 
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We have the same variational formulations, so the solution are identical. 

(ii) . Calculation of the solution X: we have Vj G {Q...N],Xa j ~ Xi j = 0. We 

reason by induction on /c. Assume that on the column A:— 1 ({(i,j); i = fc— 1}) 
and the column k ({(i, = fc}) we had X = Q. In this case, X also vanishes 
on the column fc + 1 . 
Indeed, let j £ {0, ..,N} 

(a) If j = or j ^ N then we have indeed X(^k+i,j) — because X vanishes 
on Th, the boundary of flh- 

(b) If j G {1, .., N — 1}. Let us write the equation satisfied by Xf^^ jy. 

idtX^k.^) + iAX)^k,j) = 

[dtXt^k,]) + TT? [^^{kj) - ^(fc+lj) - ^{k-l,j) ~ ^(fcj+l) - X(k,j-l)] = 

(2.15) 

Since we assumed that: X{^k_i){t) = ^(fe-i.i)(t) = 0. We reinject this 
condition in (|2.15p and we get: 

~X(k+i.j){t) = Vt. 

So we prove that if the solution vanishes on two consecutive columns, then 
the solution vanishes on all the other columns. The hypothesis being that the 
solution X vanishes on the column i = and i = 1, we indeed proved that X 
is identically null. 

(iii) . We conclude that the solution X of this problem is time continuous, (it is 

constant equal to 0) and that the jump of the solution at t = T is null, so 
c = 0. 

□ 

(vi). Discrete controllability result 



Theorem 2.22. Assume that the reference state yo defined by l|2.5p satisfies the 
non- degeneracy discrete condition (|2.19p . Let be the solution of the perturbed state 
^ andZd = yoiT). 

There exist neighborhoods V(0) C C([0,T]; Wh) and V{Zd) C J-{VLh) such that for 
all Z e V{Zd) there exists ip S V(0) such that yip{T) = Z . 

Proof. We prove this theorem using a local surjectivity result on the map Kh 
defined in l|2.7p . To this purpose, we only need to check that o?A/i(0) is surjective. 

Form theorem 12.171 dKh[Q) is surjective if: 

with 

{-dtX + AX =tr\l^j,{c) 
\x{t^Q) =0 

implies that c = 0. The proposition 12.201 shows that when the reference state satisfies 
the discrete non-degeneracy condition, then X|ri = 0. In theorem 12.211 we proved 
that this second relation implies that X = and that c = 0. We deduce that dKh{0) 
is surjective, which achieves the proof. □ 
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2.2. Semi-discrete wave equation in a square. In this section we proceed 
the same way as in the section [2?T] to prove the controllabiHty of the waves equations. 
We still consider the system in two dimensions, in the set [0, a] x [0, b] G R^, discretized 
with a step h. The domain we consider is 



n,, = {m^{ih,jh);{i,j) e 0,...,Mx 0,...,7V}. 

We already defined for m = {ih, jh) G (Zft,)^, B{m) the neighbors of m, B{m), 
the discrete interior Jl^, the boundary Th, and lastly the exterior Fh. 
We assume that the only free part of the boundary is = 0}. 

2.2.1. Unperturbed state. We define as in the parabolic case the unperturbed 
state, introducing the finite differences operator A with Dirichlet boundary conditions. 



such that 



Vm [A(t>\„^ = — 



(2.16) 



The reference state u is solution of the equation 



d^Uo + Auq 
uo(i = 0) 
[dtua{t = 0) 



= F 

= uo 

= Ul 



(2.17) 



The function F is defined the same way as the in the parabolic discrete state. 

2.2.2. The perturbed state. As in the continuous state, we are interested in 
small perturbations of the shape of the domain flh ■ We are quite free in the choice of 
the admissible transformations, and only look for sufficient conditions for the exact 
controllability. The first assumption we make on the perturbation is that the shape of 
the domain will be modified only moving nodes of the mesh using transformations 
in time. Moreover, the only moving nodes are located on the line a; = 0, and will 
move along the normal to the boundary. Because of the finite propagation speed of 
information in the wave equation, the problem will be well posed if the boundary does 
not moves faster than the information, i.e. the differential of the deformation should 
not have a module greater than the information propagation speed, in our case 1. 

On those perturbed open sets, the operator approximating the Dirichlet Laplacian 
is identical to the one defined for the heat equation. 

Definition 2.23. We consider {VfJ = l..iV - 1} the vector field fth ^ 
defined by : Vj G {I..N - 1}, 

. V,(m) = (0,0) ifm^{0,jh) 
• V,{m) = (1,0) ifm^{0,jh) 
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We denote Wh the real vector space spanned by (Vj)jg{i.. jv-i} • The admissible trans- 
formations we consider lie in the space: 

JV-l 



and such that sup ||Aj||oo < 1/2, ||9tAj||oo < 1 f 

7"=l..Af-l 



Definition 2.24. Let M{lp) denote the 2n x 2n matrix: 

-id 
Aiif) 

Definition 2.25. The perturbed state, denoted u^{x,t) e J^(ri/i), is the unique 
solution in C{[0,T], T{flh)) of the semi-discrete problem: 

^v(i-O) ={^^ (2.18) 

Proposition 2.26. We define := ^ ^ For all (f e V, U{(p) is well 

defined and lies in C([0, T], ^^($1^))^ 

Proof. [Existence et uniqueness of the solution U^p (t) :] We can indeed write the 
problem as: 




which is a Cauchy problem of the form dt~x = X{~x, t), with X{'x, t) = M{ip{t))'x + 

F{t), in R(^^~^'" « T{nh)^- The function X{~x,t) is measurable in t; it is also 
Lipschitz in Ie since A is bounded (See prop. 12. Sp . We can deduce the local existence 
and uniqueness of the solution, which in its definition domain will be continuous. □ 
Proposition 2.27. The solutions of the later problem are defined for all t > 
Proof. Gronwall's lemmas applies (See the proof of prop. I2.10p . □ 

2.2.3. Controllability of the semi-discrete wave equation. Like in the 
continuous case, we consider the function: 

A.:/^ -^^^^^ (2.19) 

where U^{T) is solution of the equation l|2.18p . Let Zd = uq{T), where uq is the 
reference state defined in l|2.17p . The question we address is the following: 

Is there neighborhoods V(0) <eV of the reference domain and V{Zd) G J-i^h) of the 
trace at t = T of the reference state such that V{Zd) C A;i(V(0)). 
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In this case we will also apply the local surjectivity theorem to A/j. 

(i) . Differentiability 

We denote here again tr the trace function a,t t ~T. The map A/j is the composi- 
tion of the function S : — > and the trace function. The Frechet-differentiability 
of A^ is equivalent to the Frechet-differentiability of 5 at 0. 

This differentiability is an immediate consequence of the differentiability of A{(p) . 
Indeed, by Cauchy-Lipschitz' theorem with parameters, if A is Frechet-differentiable 
in (fi, then the matrix 

-<-'^(.4M f ) 

is differentiable in (p which gives the Frechet-differentiability of Uip in ip. 

Proposition 2.28. cj) ^ is differentiable at in the direction ip, and the 
Prechet- differential {dAh{0),Tp), denoted Y^, is solution of the differential equation: 

\ Y^{t = Q) ^ 

(ii) . Adjoint state technique 

We have A^, = trt^T ° S. Moreover, recall that 




C{[{),T],Wh) T{dh) 



and Y^ is solution of 

r dtY^ + M{0)y^ = -{M{0),p)Uo 

1 y^(t = o) = 

Eventually, we denote L : ip — > Y^. 
We have dA,i(0) = tr\t=T ° L. 

First of all let us prove that the differential of A^ at is surjective, using the 
adjoint state method. 

dA,,(0) is surjective {c e T{dhf;'i^ G V(dAh(0)(^, c) = 0} = {0} 
On an other hand, 

Vv?eV (dA,,(G)^,c) =0 
^ V(^eV {tr\t=T{S{^)),c) ^0 

eV{S{ip)M*t=T{c)) =Q (2.20) 
Furthermore, we have S{ip:) = is solution of the differential equation: 

dtY^ + Af (0)2/^ = -Af' yo 
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Definition 2.29. The adjoint state and tr\^^j,{c) is the unique solution X 
of the equations: 

-dtX + M*X = tr|U(c) 

X{t = 0) = ^ ' 

With this definition, we replace l|2.20p tr\l^rp(c) by its expression in function of 
the adjoint state X |(23T1) . 
We deduce the following theorem: 

Theorem 2.30. The differential dAfi{Q) of hh at = is surjective if and only 
if we have the following uniqueness property: 

If c € J^i^lh) is such that 

<^'K2/"),.(o,,^,(,;,-0, V^eV (2.22) 



with 

-dtX + M*{0)X = tr\;^j.{c) 
X{t^O) = 

Then necessarily c = 0. 

(iii). Calculation of the differential of 71/ at 

We have: 

Proposition 2.31. 





(2.23) 



Af'(O) 



A'{0) 



Where A'{0) is defined in [2l8[ 
(iv). Condition X\-pi =0 



Proposition 2.32. Assume that the reference state (without the time differ- 
satisfies the discrete non- degeneracy condition on the whole boundary offth- 
Then the relation l|2.22p implies that 

X\ri = 0. 



Proof The time differential vanishes directly, and we are lead to the proposition 
[2:201 □ 

(v). Unique discrete continuation 

The aim of this section is to show the uniqueness property (|2.22p .From proposition 
12.321 under the discrete non-degeneracy condition on the reference state yo, we have 
X|ri = 0. We now show that this condition, together with the definition of the adjoint 
state l(239l) : 

-dtX + M*{0)X = tr|*^j,(c) 
X{t = 0) = 

implies that X is identically vanishing, so tr^^j, ~ and c = 0. 

To this purpose, we will study the zeros propagation of X on Qh from its boundaries. 
Theorem 2.33. The relations: 

-dtX + M*{Q)X = tr\*^r^{c) 

X\ri = 
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implies that X = and c — 0. 
Proof. 

(i) . The equation —dfX + AI{0)*X = tr\^^rp(c) can be interpreted the same way 

as the one of the heat equation. 

(ii) . The same method, returning to the equation with second derivatives, gives 

us dfu^ — which impUes, according to the initial conditions, that u^p = 
dtu^ = 

(iii) . Hence the solution X of this problem is continuous in time (it is constant 

equal to 0), so the jump of the solution at i = T is null, which imposes that 
c = 0. 

□ 

(vi). Discrete controllability result 

Theorem 2.34. Assume that the reference state ya defined in (|2.2p satisfies 
the discrete non degeneracy condition \2.1(A Let he the solution of the perturbed 
equation \2.^ . Let finally Zd = Yq{T). Then there exist neighborhoods V(0) C 

C{[Q,T]]Wh) and V{Zd) C Tidh)^ such that VZ e V(Zd)3(^ G V(0) such that 
Y^{T) = Z. 

Conclusion. In this paper we proved that the linearized heat equation was ap- 
proximately controllable with respect to the shape of the domain, while the wave 
equation is locally exactly controllable. We addressed the same questions in the case 
of the semi-discrete equations in two dimensions in a square and we proved that the 
two types of equations are exactly controllable. Nevertheless, the methods we devel- 
oped in this paper do not allow us to see the discrete control as an approximation of 
the continuous control in the wave equation. Another discretization method should 
be used to address this question, the mixed finite elements method. Indeed, we claim 
that one of the main obstacle to this interesting issue is the discretization method 
used, which does not behaves smoothly in the limit ft- — > 0. For instance we know 
(see |151 [9]) that in boundary control problem of the unidimensional wave equation, 
spurious modes with high frequency numerical oscillations appear and the observabil- 
ity constant tends to infinity when h tends to 0. It has been proved also that this 
semi-discrete model is not uniformly controllable in the limit /i — > 0. 

Nevertheless, the results of Castro and Micu in [3] are promising. They studied 
a system based on a mixed finite element space semi-discretization the linear 1-D 
wave equation with a boundary control at one extreme. They show that the con- 
trols obtained with these semi-discrete systems can be chosen uniformly bounded in 
L^(0, r) and in such a way that they converge to the HUM control of the continuous 
wave equation, i.e. the minimal L^-norm control. This result motivates to study 
in contrast to the classical finite element semi-discretization a mixed finite element 
scheme. 
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